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Abstract: We investigate the chemical (hydroxyl and ozone) and dynamical response to changing from 
present day to pre-industrial conditions in the GISS GCM. We identify three main improvements not 
included by many other works. Firstly, our model includes interactive cloud calculations. Secondly we 
reduce sulfate aerosol which impacts NOx partitioning hence Ox distributions. Thirdly we reduce sea 
surface temperatures and increase ocean ice coverage which impact water vapor and ground albedo 
respectively. Changing the ocean data (hence water vapor and ozone) produces a potentially important 
feedback between the Hadley circulation and convective cloud cover. Our present day run (run 1, control 
run) global mean OH value was 9.8x1 (f molecules/cm\ For our best estimate of pre-industrial 
conditions run (run 2) which featured modified chemical emissions, sulfate aerosol and sea surface 
temperatures / ocean ice, this value changed to 10.2xl(f molecules/cm^ . Reducing only the chemical 
emissions to pre-industrial levels in run I (run 3) resulted in this value increasing to 10.6xl(f 
molecules/cm\ Reducing the sulfate in run 3 to pre-industrial levels (run 4) resulted in a small increase in 
global mean OH (10. 7xl(f molecules/cm^). Changing the ocean data in run 4 to pre-industrial levels 
(run 5) led to a reduction in this value to 10.3xl(f molecules/cm\ Mean tropospheric ozone burdens were 
262, 181, 180, 180 and 182 Tgfor runs 1-5 respectively. 


1. Introduction 
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Identifying trends in tropospheric ozone (O3) occurring since the pre-industrial (PI) 
era is made difficult due to insufficient spatial and temporal observations, a relatively short 
lifetime (of the order weeks in the planetary boundary layer (PBL) to months in the free 
troposphere), and the fact that this specie is not emitted directly into die atmosphere. 
Therefore modeling studies are currently the preferred tool for tropospheric trend 
assessment. In recent years 3-dimensional simulations with chemistry have emerged. Of 
these, estimates of mean ozone increases range from 40% [Roelofs et al, 1997] to 60-70% 
[Wang and Jacob, 1998; Lelieveld and Van Dorland, 1995]. Changes in hydroxyl (OH) are 
more uncertain, being subject to the opposing effects of increasing CO and CH4 (which 
directly lower OH) and increasing O, (=0^P+0'D+03) and NO* (which increase OH via j03^ 
O D+H20^20H and H02+N0"^0H+N02 respectively). Estimates are in the range: 17% 
decrease [Brasseur et al.. 1998] to a 6% increase [Bemsten et al, 1997] in mean OH 
concentration since the PI era. These, and additional studies also report radiative changes due 
to O, increases. Moving from clear sky to cloudy conditions tends to lower the global mean 
forcing value by 20-60% [Haywood et al., 1998; Forster et al., 1996]. However, the above, 
and other GCM studies with chemistry [Chalita et al, 1996; Bemsten et al, 1997; Van 
Dorland et al, 1997; Levy et al, 1997] adopted climatological cloud cover when calculating 
tropospheric increase. In this study we include an online, prognostic calculation of clouds. 
We also investigate the chemical effect of reducing the sulfate aerosol distribution. Finally, 
we investigate the effect of changing our fixed monthly climatological mean sea surface 
temperatures (SSTs) and ocean ice distribution from present day (PD) to PI . 

We employ a 3D model, the Goddard Institute for Space Studies (GISS) general 
circulation model (GCM) in order to gain new insight into anthropogenically-driven oxidative 
and radiative changes to the troposphere from PD to PI conditions. In summary, compared 
with the PD control run, we have modified to PI conditions: (1) chemical emissions, 
including the ozone precursors (2) the sulfate aerosol distribution (3) the oceanographic data 
and (4) the greenhouse gases (GHGs). We report resulting changes to the tropospheric 
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hydroxyl (OH), O, (= 03 + 0 'D+ 0 ^P) and water vapor concentrations as well as the associated 
dynamical/cloud response. Section 2 provides a model description and overview of the runs; 
section 3 discusses results; section 4 presents conclusions and ftiture work. 


2. Model description 

The GISS n’ (two prime) GCM [Hansen et al., 1997] was developed from the GISS 
GCM version II which is described fully by Hansen et al. [1983]. GISS II’ is a gridpoint 
model having 72 and 46 longitudes and latitudes respectively (-4x5 degrees horizontal 
resolution) and 9 sigma coordinate levels in the vertical with a dynamical top at lOmb. The 
vertical resolution is admitedly crude which has repercussions both for modeling ozone and 
the prognostic cloud calculations. This point is discussed further in section 4 (conclusions). 
For the PD run, fixed climatological monthly mean Sea Surface Temperature (SST) data 
based on the Atmospheric Modeling Intercomparison Project (AMIP) 1979-1993 and ocean 
ice coverage are incorporated. Tracer advection is achieved via the quadratic upstream 
method [Prather, 1996] which produces little numerical noise and realistically maintains 
strong tracer gradients. Further details may be found in the companion paper in Part I, 
Shindell et al., hereafter referred to as Part I. 

We employ 10 chemical tracers (O*, NO*, N 20 j, HNO 3 , H 2 O 2 , CH 3 OOH, HCHO, 
HO 2 NO 2 , CO and H 2 O) using the ‘family’ approach where appropriate with 37 reactions for 
24 molecules [DeMore et al, 1997]. Peroxyacetyl nitrate (PAN) and explicit NMHC 
chemistry is omitted. A proxy for CO (60 Tg C year * i.e. 10% of the isoprene burden) based 
on [Wang et al., 1998] is included. Other NMHC chemistry is omitted. We therefore tend to 
underestimate ozone production, especially in the PEL. Our ‘end product’ proxy approach 
for isoprene ignores ozone production via the intermediate oxidation channel RO 2 +NO. 
Other than isoprene, ‘equivalent CO’ for other NMHCs was not included. The FASTI 
photolysis scheme, after validation with results from the Photochemical Model 
Intercomparison (PhotoComp); [Olson et al, 1997], is employed. The chemical and 
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dynamical timesteps are 1 hour, 5 minuter respectively. N 2 O 5 hydrolysis on sulfate aerosol is 
calculated according to Dentener and Crutzen^ [1993] using monthly-mean modeled sulfate 
surface area taken from Koch et al [1999]. Wet and dry deposition are adapted from Wang et 
al [1998]. The lightning parametrization generates NO^ interactively with the convection 
scheme [Price et al., 1997]. Other chemical emissions are prescribed as monthly-mean values 
and include CO (biomass, equivalent isoprene and industrial) and NO, (industrial) adapted 
from Wang et al., [1998]; NO, aircrafr fleet emissions for 1992 [Baughcum et al., 1996; 
Schulte et al., 1997] and NO, soil emissions adapted from the GEIA l^xl" database. Emission 
inventories are climatological where possible. Further details may be found in Part I. As 
described there, the model does a good job of reproducing PD observations of 0„ and based 
on the lifetime of methyl chloroform, yields a good global OH field. 

All runs begin December 1 ** and were ten years in duration. Initial analyses suggested a 
spin-up time of 1 year and 2 years for the control run and for all remaining runs respectively. 
Data obtained during spin-up were subsequently discarded. All figures show January and July 
zonal mean values unless otherwise indicated. Table 1 gives an overview of the runs. Run 1 is 
the PD control run. Run 2 is our best estimate of PI conditions. It featured chemical 
emissions, GHGs (CO 2 , CH 4 and N 2 O) and sulfate aerosol reduced from PD to PI conditions. 
Industrial and aircraft CO and NO* emissions were set to zero; methane in the chemistry 
scheme was reduced from 1.803 ppbv (SH=1.750, NH=1.855) to 0.730 ppbv in both 
hemispheres; biomass burning emissions were set to 10% of their PD value but retaining the 
same geographical distribution [Wang and Jacob, 1998]and sulfate aerosol was changed to a 
PI distribution based on the GCM sulfur-cycle model described in Koch et al. [1999]. These 
implied veiy large, widespread increases (up to 5000%) in sulfate since the PI era over 
industrialized areas, especially in NH summer. Run 2 also featured specified climatological 
monthly mean SSTs and ocean ice data but these were modified from PD to PI (1850) 
conditions. The 1850 ocean data were obtained from a previous integration of the GCM at 
8 x 10 “ horizontal resolution which calculated rates of ocean heat fluxes (the ‘qflux’ scheme) 
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for the PI atmosphere [Hansen et al, 1983; Russell et ah, 1985]. Part 1 gives further 
information. 

Runs 3-5 are sensitivity runs (see Table 1). Run 3 most closely resembles other works 
in the literature investigating ozone change since the PI era in that it features PI chemical 
emissions but has PD sulfate aerosol and ocean data. Run 4 tests the effect of reducing sulfate 
aerosol in run 3 to PI conditions. Run 5 tests the effect of changing the ocean data in run 4 
from PD to PI conditions. When calculating O, burdens in the troposphere we defined the 
tropopause using a cutoff criteria of 150ppb ozone. Changing this value by +/- 10% led to 
changes in the burden of +3/-4% respectively. 

3. Results and Discussion 

Table 2 shows a complete overview of the global tracer chemical budget. The values in 
brackets represent (PD-PI)/PI expressed as a percentage. O* (obtained separately as the 
tropospheric burden) and global OH concentrations are discussed separately in the sections 
which follow. Firstly, a general feature in Table 2 is enhanced wet deposition for the warmer, 
damper PD atmosphere. Secondly, all tracers show increases in the PD run occurring mainly 
in the NH. The increases in NO^, HNO3, HO2NO2 and CO are almost entirely associated with 
changes in precursor emissions (run 3). Increases in H2O2, CH3OOH and HCHO since PI are 
also chiefly associated with precursor emission changes but also have an added component 
(-5-15%) associated with warming the ocean/moistening the atmosphere (run 5). 

Ox chemical production is negative in all hemispheres for all PI runs but strongly 
positive for the PD run in the NH. The cold ocean (run 5), with 13% less water vapor than 
the PD run, features around 20% more Ox chemical production than other PI runs not having 
the cold ocean, due to a weaker 0 'D+H 20 sink. The lowered sulfate run (run 4) led to a 
significant slowing in the global mean rate of N 2 O 5 chemical destruction, from - 0.8 to - 0.66 
Kg/s in runs 3 and 4 but otherwise had little impact. However, more significant effects are 
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apparant when zonal mean difference plots are obtained (section 3.2). Global lightning 
increased by 8% in going from the best PI (run 2) to the PD, 7% of which was associated with 
changing the ocean data (run 5). Global water vapor similarly increased by 14%, 13% of 
which was associated with changing the ocean data. Changes due to chemistry became 
generally stronger (i.e. stronger destruction of NO, and CO, stronger production of die other 
tracers) in going from PI to PD. The exception was HOjNOj which featured overall weaker 
chemical destruction for the PD run (e.g. by 52% compared with the best PI, run 2), and 
displayed increased PD production in the NH associated with increased NO, yet increased 
destruction in the SH proabably associated with enhanced thermal decompostion in the 
warmer PD atmosphere. Changing to PI ocean data had only a small effect (a couple of 
percent at most) upon the chemistry terms in run 1, except for 0„ which increased by 
around 20% due to the suppressed 0 'D+H 20 sink, as already mentioned. 

Table 3 asseses model internal variability by presenting confidence intervals at 99%, 
95% and 60% calculated from the t-distribution based on eight Julys for NO,, NO, from 
lightning, OH, H 2 O and tropospheric ozone. Data is shown for run 4 although all PI runs 
displayed similar internal variability. Table 3 data are considered in the sections which follow, 
in which we mainly discuss changes due to ozone and OH in the zonal mean for the various 
runs. 

An interesting application of coupled chemistry-climate models is to quantify 
changes in the downward fluxes of O, and NO, from the stratosphere to the troposphere. 

Table 4 presents vertical mass fluxes for these two families between level 7 and level 8 (mean 
pressure=110 mb). TTie same definition was applied in Part 1. The values in brackets 
represent (PD-PI)/PI expressed as a percentage, as for Table 2. Results in Table 4 suggest 
that in going from PI to PD, the net downward transport of O, from the stratosphere has 
generally slowed (e.g. by 3.6% for the best PI estimate, run 2) whereas that for NO, has 
increased (e.g. by 37.5% for run 2). Further, the O, changes appear to be driven mainly by 
changes in the emissions (run 3) and occur more strongly in the NH whereas the NO, changes 
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are driven mainly by changing the ocean data (run 5) and occur more strongly in the SH. 
Since NO, and O* are fixed and proportional to each other in the stratosphere, the differing 
response reflects changes in the upward, rather than the downward fluxes. A slowing in net 
downward O* in Table 4, controlled mainly by emissions, is consistant with higher values of 
0„ hence an increased upward flux in PD. An increase in net downward NO*, controlled 
mainly by the ocean data, is consistant with a damper atmosphere leading to faster loss of 
NO* into HNOj in PD, hence a slower upward flux. Changing the sulfate aerosol has only a 
small effect (compare runs 3 and 4). 

3.1 Control (run 1) versus Best PI (run 2) 

Figure 1 shows O* (ppbv) for the control (run 1) while Figure 2 shows the difference 
in O* (ppbv) for run 1 and run 2 (best PI) (annual averages were shown in Part I). The total 
tropospheric O* burden for runs 1 and 2 (Table 5) was 262 and 181 Teragrammes (Tg) 
respectively i.e. an increase of 45% (significant to > 99%, Table 3) fi’om PI to PD. Peak 
differences of 30-40 ppbv occur in the NH during July throughout much of the troposphere. 

The winter hemisphere in Figiue 2 displays anomalous, negative (5-10 ppbv) (PD-PI) 
O* features near the tropopause. Figure 3 shows one standard deviation O* of eight Julys for 
the control run on levels one and seven. January values (not shown) were generally lower by 
around 10% in the NH and were similar in the SH. Standard deviation values on the model 
sigma levels are enhanced over high terrain e.g. Greenland and Antarctica where absolute O* 
values are high. Also, level seven (202mb) values are much higher (20-50 ppbv) than those 
on level 1 (959mb) which are typically 0-15 ppbv. The negative O* features near the 
tropopause are only of the order lo. However, note that the radiative forcing is most 
sensitive in regions where these features occur and there does exist a possible physical 
explanation for their occurrence: NO* differences (ppt) for (control PD and best PI) i.e. (run 
1 and run 2) (not shown) indicated higher PI NO* near the tropopause but not exactly at 
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those latitudes corresponding to the Ox anomalies. Further investigation discounted 
stratospheric-tropospheric exchange as being the source of the high NO* features. The 
chemistry scheme was the main contributor. The features corresponded reasonably well (but 
not exactly) to PI HNO3, which was lowered in these regions by roughly an equal magnitude. 
This then suggests that lowered PI OH was playing a role. The features were not exactly co- 
located since HNO3 has a longer lifetime than NO, and is affected by lar^e-scale dynamics to 
a greater extent. The sign of (PD-PI) OH change in the real atmosphere is uncertain due to 
opposing chemical effects, as discussed in the introduction. Lightning in the tropics also 
played a lesser role (not shown) occasionally leading to higher PI NO, in some regions but 
more usually resulting in lowered NO*. Run 1 (control) and run 2 (best PI estimate) featured 
3.87 (NH=2.33, SH=1.55) and 3.57 (NH=2.08, SH=I.49) Tg N / year respectively. Such a 
lowering of 8% is significant at the 99% confidence interval (Table 3). The bulk of the 
decrease was due to changing the ocean data from PD to PI conditions, as will be shown in 
section 3.2.3. Another possible cause of high PI ozone is enhanced clouds leading to a 
lowered jO'D sink. However, clouds only increased in a narrow latitude belt centred over the 
equator (Figure 13); since the high PI ozone occurred mainly at mid to high latitudes clouds 
were concluded not to play a role. 

Figure 4 shows OH (1x10* molecules/cm*) for run 1 while Figure 5 shows the 
difference (1x10* molecules/cm*) for OH (run 1 - run 2) i.e. (PD control - best PI). Annual 
averages were shown in Part I. Global mean OH values for runs 1 and 2 (Table 3) were 9.8 
(NH=12.0, SH=7.6) x 10* molecules/cm* (almost the very same value of 9.7 was derived by 
Prim et al, [1995] using CH3CCI3 measurements between 1978 and 1994) and 10.2 
(NH=12.0, SH=8.4) x 10* molecules/cm* respectively i.e. a decrease of 3.9% (NH=0.0%, 
SH=9.5%) from PI to PD, which is significant at the 99% confidence interval (Table 3). PI 
OH is generally higher than PD throughout the tropical troposphere except at (1) higher 
levels especially in summer and (2) near the ground especially between the tropics and mid- 
latitudes. Higher PI OH is associated with lower CO and occurs despite lower NO, and lower 
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water vapor; PD water column (mm) = 22.51 (NH=23.11, SH=21.92) whereas PI water 
column = 19.49 (NH=19.71, SH=19.27). 

Figtire 6 shows the difference (tenths "C) in temperature for (PD control - best PI) 
i.e. (run 1 - run 2). Stratospheric cooling leads to a 1-2°C lowering in temperature; tropical 
mid-troposphere heating results in 3-5°C warming. Annual mean surface tem]>eratures 
increased from ILS^O^C (NH=12.546, SH=10.234) in run 2 to 13.343®C (NH=14.333, 
SH=12.354) in run 1, corresponding to 17.1% (NH=14.2%, SH=20.7%). Note that for the 
control run, annual mean sigma = 0.025K; this value is somewhat lower than observations 
which are in the range 0.2-0.4K [IPCC, 1996]. The discrepancy is associated with fixing SSTs 
in the model. 

Global mean air temperature increases were 1.853°C (NH=1.676, SH = 2.030). These 
were associated mainly with changing the ocean data and were stronger in the SH due to a 
larger area of ocean coverage. This more than outweighed greater O, changes in the NH. 

The difference (%) in moist convective cloud cover (not shown, but very similar to 
Figure 1 1) for run 1 (PD control) - run 2 (best PI) revealed that the equatorial mid 
troposphere is more cloudy by 2-4% in the PD run (aj„=0.50, ajui=0.53% at 323mb). This 
region is straddled north and south by features for which the opposite is true. We interpret 
this results as follows: increased PD water vapor in the central fropics implies extra radiative 
heating hence decreased stability which strengthens the ascending Hadley cell hence increases 
moist convective cloud cover. The adjacent, descending arms are also strengthened, which 
tends to decease cloud cover in these regions. The changes arise mainly due to changes in the 
ocean data and also impact upon NO* from the lightning scheme. They are quantified and 
discussed frirther in section 3.2.3, which discusses the effect of the ocean data alone. 

3.2 Effect of precursor emissions, sulfate aerosol, ocean data and GHGs 



11 


The only modification to run 3 compared with the PD control (nm 1) is a lowering in 
O, precursor emissions. Run 3 is therefore the most comparable with other works which 
simulate O, increases since the PI era, as discussed in 3.2.1. In 3.2.2 we investigate the effect 
in run 3 of modifying to PI conditions the sulfate aerosol in the chemistry scheme (run 4); in 
3.2.3 we then investigate the effect in run 4 of modifying to PI conditions the ocean data 
(run 5). The best-estimate PI run (run 2) is then obtained by modifying to PI conditions the 
GHGs in run 5. In 3.2.4 we compare run 2 (best PI) with run 5 which differ only in that run 2 
has PI levels of GHGs. 

3.2.1 Precursor emissions. Figure 7 shows the difference (ppby) for O, run 1 (PD control) 

- run 3 (PI emissions i.e. most comparable to other works). Total O, burdens (Tg) for run 3 
(Table 5) were 180 Tg, an increase of 45% from PI to PD. Our % change value compares 
reasonably well with other works e.g. (360.0 to 230.4 Tg) i.e. 56.3% {Mickley et al. 1999]; 
(273.6 to 192.0 Tg) i.e. 42.5% [Roelofs et al. 1997] and (312.0 to 192.0 Tg) i.e. 62.5% 
{Wang and Jacob, 1998] (see also discussion in Part I). In Table 5, the ozone change for run 
3 (PI emissions) is similar to that of run 2 (best PI) and so the precursor emissions are 
responsible for the bulk of the ozone change. Figure 7 featured an anomalously high PI O, 
feature occurring in SH high latitudes (but less prevalent than for the best PI run, Figure 2). 
This arose due to higher PI NO^ (as for run 2, section 3.1) and was associated with lower PI 
OH (not shown); again its magnitude was close to one standard deviation in O, (Figure 3). 

Plotting the difference (IxlO’ molecules/cm^) for OH run 1 (PD control) - run 3 (PI 
emissions) revealed results (not shown) similar to Figure 5, which showed run 1 (PD control) 

- run 2 (best PI estimate). The global mean OH value for Run 3 (Table 5) was 10.6x10* 
molecules/cm* (NH=12.5, SH=8.8) corresponding to a 7.5% (NH=4.0, SH=13.6) decrease 
from PI to PD, which is significant at the 99% confidence interval (Table 3). Values from 
other works vary between a 20% increase to a 6-7% decrease (see e.g. Wang and Jacob. 1998 
or Part I). Wang and Jacob [1998] calculate similar OH changes yet do not calculate the 
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anomalous high PI NO, and O* features occurring in SH high latitudes (eis discussed in 3.1). 
Hauglustaine et al. [1994] calculate widely different OH spatial changes and also do not 
calculate the higher PI O,. Other works report only O*, not OH [Roelofs et al, 1997; 
Stevenson et al, 1998; Brasseur et al, 1998; Bemsten et al, 1997] and they too do not 
report higher PI Ox in SH high latitudes. Mickley et al [1999] do calculate higher PI NO, in 
this region (L. Mickley, personal communication). However they do not see higher PI O, 
values and have suggested that the NO, phenomenon is related to PAN formation (not 
present in our model) being suppressed in remote regions. Grewe et al [2001a] compared 
GISS modeled NO, near the tropopause with observed data and found that modeled NO, 
reasonably reproduced observations for zonal mean values, but the modeled values failed to 
capture observed longitudinal variability, possibly reflecting inadequacies in modeled nitrogen 
oxide chemistry and /or dynamical effects. 

3.2.2 Effect of Sulfate aerosol. Plotting the difference (ppbv) in O, for run 4 (PI emissions 
and PI sulfate) - run 3 (PI emissions only), (not shown) revealed mostly small values < 2 ppb 
with higher (up to 10 ppb) values occiuring near the tropopause, although these were not 
statistically significant due to higher concentrations of ozone in this region, as discussed 
previously. Negative O, features occurred in the SH in both winter and summer. Overall, 
ozone was lowered by 2-3% in the upper troposphere. Figure 8 shows the difference (ppt) for 
NO, run 4 (PI emissions and PI sulfate) - run 3 (PI emissions only). Apparent are the 
succession of positive/negative features near the tropopause. Longer runs are required to 
establish whether these features have a physical basis or arise via internal model variability. 
Total O, burdens for run 4 (Table 5) were 180 Tg, similar to run 3. 

Figure 9 shows the difference (IxlO’ molecules/cm^) for OH run 4 (PI emissions and 
PI sulfate) - run 3 (PI emissions only). Although global mean OH values were very similar for 
both runs (Table 5), values did vary quite considerably latitudinally. The complex behavior in 
the NH we interpret as sulfate on the one hand lowering OH (lowered NO, suppresses the 
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reaction H 02 +N 0 '^ 0 H+N 02 ) and on the other hand increasing OH (increased HNO 3 
photolysis). The vertical band structure is believed to arise due to changes in the Hadley cell 
circulation which results in alternately enhanced/depleted regions of cloudiness, which results 
in decreased/increased UV below cloud height as well as producing modifications in OH 
concentrations. 

3,2.3 Effect of ocean data. 


Plotting die difference (ppbv) for 0» run 5 (PI emissions, sulfate and ocean) - run 4 
(PI emissions and sulfate only) revealed a positive feature (up to 20 ppb) occurring in the 
tropical upper troposphere during July (not shown). Similar to the features discussed in 3.2.1, 
its magnitude is only of the order one standard deviation (Figure 3) so it may have arisen via 
model variability. However, there is a possible physical explanation, namely higher NO, (not 
shown) in run 5. This was associated with lower OH (3.9% i.e. significant at the 99% 
confidence interval, Table 3) in run 5 (Figure 10), which was in turn associated with lower 
water levels (not shown) as a result of lower evaporative fluxes fi'om the colder PI ocean. 

Total tropospheric O, burdens for run 5 (PI emissions, sulfate and ocean) in Table 5 
were 182 Tg corresponding to a 1 . 1 % increase (of the order 95% significant in Table 3) on 
cooling the ocean. This arose despite a concomitant lowering in NO, input from the 
lightning scheme fi'om 3.86 (NH=2.32, SH=1.53) TgN/year to 3.62 (NH=2.12, SH=1.51). As 
already noted, our PI ocean run features lower HjO hence the 0 'D+H 20 sink is slowed which 
may account for the increase in PI Ox, at least in part. Consistent with this result, 

Fuglestvedt et al. [1995], Johnson et al. [1999] and Brasseur et al [1998] predict a future 
decrease in ozone associated with a damper atmosphere leading to a faster 0 'D+H 20 sink. 
[Grewe et al, 2001b] however, employed future (2015) SST data in the GISS model and found 
increased NO*, hence ozone (of the order 10 % in the tropics) which they attributed to 
changes in circulation, hence precipitation and washout. For the runs here, however, NO, 
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changes only very slightly when we cool the oceans (Table 2). Global mean OH values for run 
5 were 10.3x10* molecules/cm* (NH=12.2, SH=8.5) corresponding to a 3.9% (NH=2.5, 
SH=3.5) increase in changing from PI to PD ocean data (run 5 compared with run 4), which 
is significant to 99% (Table 3). 

On plotting the change (%) in moist convective cloud cover due to cooling the ocean 
(Figure 11), run 5 (PI emissions, sulfate and ocean) - run 4 (PI emissions and sulfate only) it 
was clear that the bulk of the cloud behavior noted in the best PI (section 3.1) run was 
attributable to changes in the ocean data. We now quantify this effect. Lower moisture fluxes 
from the colder ocean leads to decreased tropical water vapor (~20%), decreased radiative 
heating hence decreased (10-16 x 10* jW/sigma layer) available potential energy (PE), 
slower ascent of the Hadley cell and decreased convective cloud cover (e.g. by 2-4% in July in 
the tropics, 0=0.53%). Slowing the ascending Hadley cell also slows the adjacent, descending 
cells of dryer air, which leads to an increase in convective cloud cover in these regions, also 
in the range 2-4%. 

Large increases in available PE (30-40 J/mVsigma layer ) were also apparent over 
Antarctica diuing PI, possibly due to an enhanced ocean ice-albedo effect, although no 
corresponding increase in cloudiness was observed. This requires further investigation. Figure 
12 shows changes in the mass streamfunction, (MS), for run 5 (PI emissions, sulfate and 
ocean) - run 4 (PI emissions and sulfate only). The MS is an indicator of the mean meridional 
circulation. Positive (negative) values indicate anti-clockwise (clockwise) motion in the plane 
of the Figure. In Figure 12, zonal mean rate of vertical overturn (lO’ kg/s) for the control 
run, MSOgiob,! = 1.40, 0.97x10® kg/s in January, July respectively while at 550mb, 14S MS a = 
10.28, 6.47x10® kg/s for Jamxary, July respectively. Figure 12 reveals decreased upward 
motion in the tropics and reduced downward transport of dry air in the subtropics, consistent 
with the cloud response. Corresponding features were noted in the input of NO* by the 
interactive cloud-lightning scheme. 
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Inspecting the difference (®C) in temperature for run 5 (PI emissions, sulfate and 
ocean) - run 4 (PI emissions and sulfate only) (not shown) revealed that changing the ocean 
data accounted for the bulk of the temperature change in the best estimate PI run, i.e. results 
were almost identical to Figure 6, which showed temperature for control - run 2 (best PI). 
The drier PI troposphere is colder by 2-4°C over widespread regions, especially in the tropics 
aroimd 300mb. The drier PI stratosphere warmed by around 1®C. Global, mean temperature 
(a=0.025°C) was higher by 1.629“C (NH=1.359, SH=1.900) for run 4 (PD ocean) compared 
with run 5 (PI ocean). (Here we quote global mean temperature to three decimal places 
because this order of accuracy is suggested by the standard deviation of model variability; we 
do not imply that die model reproduces this level of accuracy compared with observations). 

3.2.4 Effect of GHGs. 

On plotting the difference (ppbv) in O* for run 2 (best PI estimate) - run 5 (all PI 
changes except GHGs), again we find small, anomalous higher PI O, features (not shown) 
near the tropopause which are close to one standard deviation (Figure 3). Also, there exists a 
seasonal pattern over the pole at around 200mb with lower (higher) O, in the summer 
(winter) for run 2 compared with run 5. Changes in NO* (not shown) mosdy explain the NH 
Ox behavior but not the SH. Figure 13 shows the difference (tenths °C) in temperature (run 2- 
run 5). There features weak (0.2-0.5 °C) cooling during PI over wide areas in the 
troposphere, which may not be statistically significant (e.g. s at 200mb in July=0.53°C). 

There is a possible, albeit weak, chemical contribution, namely the reaction HO 2 +NO which 
mediates in-situ chemical ozone production and which shows a wcjddy negative temperature 
dependence, typical of many radical recombination reactions. Thus, cooling is expected to 
lead to enhanced ozone via this process, although the weakness of the temperature 
dependance suggests ozone changes of only a couple of percent at most. 
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There also features PI cooling of up to one degree in high latitudes in the SH winter. 
Polar regions in the dark are sensitive to dynamical changes which affect polewards transport 
of heat. Lowering the GHGs led to a decrease in the rate of descent of the polar cell, 
consistant with the cooling feature. The phenomenon only features as a response to the 
GHGs being lowered; it is not evident e.g. for the best estimate PI run. 

The total tropospheric O, burden for runs 5 and 2 (Table 5) were similar, 
corresponding to 182 and 181 Tg respectively. O, mixing ratios however, tended to change 
most near the tropopause and may therefore impact the radiative calculations to a greater 
extent as will be discussed in section 3.3. Global mean OH values for run 5 and run 2 were also 
quite similar, corresponding to 10.3x10* molecules/cm* (NH=12.2, SH=8.5) and 10.2x10* 
molecules/cm* (NH=12.0, SH=8.4). 

Global mean temperature increases were 0.153°C (NH=0.231, SH=0.076) for (run 5 - 
run 2) i.e. (PD GHGs - PI GHGs); a=0.025"C. These values correspond only to reductions in 
the concentrations in air of the GHGs. They do not reflect feedbacks associated with ocean 
heating. 

3. 3 Radiative feedbacks: chemical composition and clouds 

For any given, instantaneous forcing, the resulting response varies considerably between 
models. Over sufficiently long periods, the response will feed back upon evolving chemical 
constituents such as ozone and it will eventually impact future forcing values. Therefore, 
radiative forcing (RF) definitions must also stipulate which, if any model-dependant 
response(s) are included. The most common is that defined by {IPCC, 1996], which quotes 
RF at the tropospause after the stratosphere has been allowed to respond dynamically. 
Important quantities in our model which potentially feed back upon the radiative forcing 
include dynamical, cloud and albedo effects. Also, since we fix the SSTs our model cannot 
fully respond over the longer oceanic timescales although we believe results are valid over the 
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shorter (ten year) period employed. RF values as defined by [IPCC, 1996] are unavailable for 
the runs discussed here since the present calculations included the model-dependant response 
to ozone chmges. However, future runs planned will calculate RF differences for PD and PI 
ozone distributions. 

For run 3 (PI emissions), at NH mid-latitudes (~500mb) we calculate a region of high 
(40-50ppbv) Ox (Figure 7) which lies directly above lower values. In the.SH this does not 
occur and there even is some indication of the reverse process taking place in the lower to 
mid troposphere at mid to high latitudes. Global mean RF increases due to enhanced 
tropospheric ozone occuring since PI is 0.4 (+/- 0.2) WW. This is significant when 
compared with the same quantity summed for (COj+CH^+NjO) =» 2.5 (+/- 0.3) W/m^ [IPCC, 
1996]. Also, vertical (and horizontal) distributions of ozone tend to be less homogeneous 
compared with the other GHGs, due to its shorter lifetime. Introducing a heating source (O,) 
high (low) in the troposphere such as the case in the NH (SH) tends to increase (decrease) 
atmospheric stability hence will decrease (increase) convective cloud cover (not shown) and 
so increase (decrease) SW absorption in the NH (SH). Models which employ climatological 
cloud cover will tend to overlook this phenomenon. 


3.3.2 Sulfate aerosol Sulfate aerosol as specified in the chemistry scheme may 
either increase or decrease NO* (hence O*) depending on location, as discussed in section 
3.2.2. For the runs discussed here, the radiation scheme assumes a PD distribution of sulfate 
aerosol (as well as dust sea-salt, organic carbon etc.) whose concentrations are slightly 
different fi-om that used in the chemical calculations. The impact of this difference is likely 
to be small. Note that only sulfate in the chemistry subroutines was reduced to PI levels; the 
radiation sulfate was left unchanged. Globally, O* and hence cloud responses were relatively 
small (2% decrease and 0.03-0.05% increase respectively). 



18 


3 J.3 Ocean data There featured relatively large responses in cloud cover to the PI 
ocean (Figure 11) which were different in sign between hemispheres (lower/higher in the 
NH/SH respectively). Clouds were overall sparser in the NH because the large, central region 
of lowered cloudiness in Figure 11 was displaced across the equator toward the NH. PI ocean 
data led to lowered (~20%) water vapor in the tropics associated with slower Hadley 
overturning and increased atmospheric stability. There appears to exist radiative feedback 
between changes in radiative water, ocean ice and cloudiness. Ozone values also changed near 
the tropopause but a longer run is required to establish significance. 

4. Conclusions 

An important result is the cloud feedback mechanism. To the author’s knowledge, other 
GCMs/CTMs which calculated O* and OH changes since the PI era mostly employed 
climatological clouds hence may produce different dynamical and cloud responses (which feed 
back upon the chemistry) compared with our runs which calculate cloud cover interactively. 
The contribution of clouds to the earth’s radiation budget is believed to fall generally within 
the range 20-60 WW for the LW heating component, with generally higher values 
associated with the SW cooling [Ramanathan, 1987] contribution. Relatively small changes 
in global mean cloud cover therefore may produce radiative changes which are significant 
compared with the radiative changes associated with anthropogenic O, increases which fall in 
the range 0.2-0.6 W/m^. Uncertainties associated with the cloud-radiation scheme and cloud 
feedback effects [Cess and Potter, 1987] are generally large in comparison. Changing the 
ocean data led to lowered H 2 O hence a slowing of the Hadley overturning which resulted in 
adjacent regions of depleted/enhanced tropical cloudiness in regions where upward/downward 
flow of moist/dry air are slowed. Therefore, both the vertical position and magnitude of the 
change in GHGs play a role. This mechanism is also relevant for future predictive climate 
runs in which the vertical distributions of the GHGs H 2 O and O 3 are expected to change. It is 


19 


unclear whether future increases in ozone will occur mainly in the lower or upper 
troposphere. Also, other chemistry models currently used to predict future climate conditions 
are mainly CTMs which lack our cloud feedback mechanism [IPCC, 2000]. The mechanism 
may be sensitive to the vertical resolution. Earlier experiments with the GISS GCM using a 
finer vertical resolution and a higher model lid (23 layers, extending to to 8Skm) tended to 
magnify moist convective fluxes and increased release of latent heat in the troposphere 
[Rind and Lemer, 1996; Rind et al., 1998] which resulted in additional upper level cloud 
cover. This will tend to reinforce the conclusions of the present work. Cloud-chemistry 
interaction is an important area of future study. However it is beyond the scope of the 
present paper to address this topic in full - in a forthcoming paper (Grenfell et al., in 
preparation) we include a complete analysis of cloud-ozone feedback. Compared to most 
other models, the GISS GCM possesses a generally simpler chemistry scheme yet a more 
detailed cloud-radiation treatment. Despite this, our chemistry scheme has been found to 
compare favorably with other, generally more complex schemes. This paper, and Part 1 have 
shown that our calculated PD and PI burdens compared well with other published works. 

We also investigated the effect of lowering sulfate aerosol to PI conditions. Kiehl et 
al. [1999] also report lowering the sulfate burden but it is unclear whether they did this 
uniformly or adopted a new PI distribution, as we did. Lowering sulfate in our model resulted 
in a succession of higher/lower NO, (and O,) features near the tropopause which require 
longer runs to establish whether they are significant or due to model variability. The overall 
effect on the O, budget was small. However, since we do not include the effects of NMHCs, 
in-situ ozone changes in our model are less sensitive to changes in NO, arising from changes 
in sulfate. 

On lowering only the precursor emissions, we obtained an ozone response similar to 
other works which lends confidence to our model results. Finally we investigated the effect of 
lowering the GHGs. This strengthened the meridional temperature gradient and stimulated 
polewards transport of heat, resulting in extra heating of up to 1 degree in SH polar night. 


r 
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Understanding the changing face of chemistry-climate interaction remains a challenging 
problem. 
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Figure Legends 

Figure 1: O, (=03+0'D+0^P), ppbv for the control run. Note that the top model level is at 
lOmb 

Figure 2; O,, ppbv for the control - best PI run (run 2) 

Figure 3: One standard deviation O* (ppbv) in July, levels 1 and 7 

Figure 4: OH (1x10* molecules/cm*) for the control run 

Figure 5: OH (1x10* molecules/cm*) for the control - best PI run^nin 2) 

Figure 6: Temperature (tenths °C) for the control - best PI run (run2) 

Figure 7: O*, ppbv for the control - pre-industrial emissions run (i.e. run 3, most comparable 
to other works) 

Figure 8; NO, (ppt) for run 4 (PI sulfate and emissions) - run 3 (PI emissions) 

Figure 9: OH (1x10^ molecules/cm*) for run 4 (PI sulfate and emissions) - run 3 (PI 
emissions) 

Figure 10: OH (IxlO'* molecules/cm*) for run 5 (PI sulfate, emissions and ocean) - run 4 (PI 
sulfate and emissions) 

Figure 1 1 : Moist convective cloud cover (%) for run 5 (PI sulfate, emissions and ocean) - run 
4 (PI sulfate and emissions) 

Figure 12: Mass streamfunction (IxlO’ kg/s) for run 5 (PI sulfate, emissions and ocean) - run 
4 (PI sulfate and emissions) 

Figure 13: Temperature (tenths “C) for run 2 (PI sulfate, emission, ocean and GHGs) - run 5 
(PI sulfate, emissions, ocean) 
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Table 1: Overview of runs 
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Table 2: Overview of global, annual mean chemical tracer budget 
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Run 1 (PD 

Uniu 

Process/ Cone. 

Global 

NH 

SH 

control) 






Ox 

10^ Kg/s 

Dry dep. 

-36.1 

-23.7 

-12.5 

Ox 

10^ Kg/s 

Chemistry 

1.2 

1.3 

-0.07 

NO, 

10 '^^ ppv 

Mixing ratio 

35,8 

40,9 

30.8 

NO, 

10" Kg/s 

Dry dep. 

-168.0 

-146.1 

-21.9 

NO, 

10" Kg/s 

Chemistry 

-1,046.7 

-842.2 

-204.6 

NO, 

10" Kg/s 

Lightning 

122.5 

73.6 

48.9 

NO, 

10 " Kg/s 

Surface 

1,072.7 

898.5 

174.2 

HNO 3 

10 " ppv 

Mixing ratio 

155.6 

182.2 

128.9 

HNO 3 

10 " Kg/s 

Moist convect. 

-928.2 

-564.6 

-363.6 

HNO 3 

10" Kg/s 

Wet dep. 

-898.3 

-774.6 

-123.7 

HNO 3 

10 " Kg/s 

Dry dep. 

-1,949.9 

-1,655.5 

-294.4 

HNO 3 

10" Kg/s 

Chemistry 

4,678.4 

3,757.3 

921.1 

H 2 O 2 

10 ‘" ppv 

Mixing ratio 

106.7 

119.3 

94.2 

H 2 O 2 

10" Kg/s 

Moist convect. 

-5,557.2 

-2,939.0 

-2,618.2 

H 2 O 2 

10 " Kg/s 

Wet dep. 

-1,064.7 

-579.8 

-484.9 

H 2 O 2 

10" Kg/s 

Dry dep. 

-278.1 

-156.0 

- 122.0 

H 2 O 2 

10" Kg/s 

Chemistry 

11,379.0 

6,045.0 

5,334.0 

N 2 O 5 

10‘^ Kg/s 

Chemistry 

93.0 

95.0 

- 2.0 

CH 3 OOH 

10 " ppv 

Mixing ratio 

24,9 

24.7 

25.1 

CH 3 OOH 

10" Kg/s 

Dry dep. 

-737.2 

-363.1 

-374.1 

CH 3 OOH 

10" Kg/s 

Chemistry 

757.8 

337.8 

420.0 

HCHO 

10 '*^ ppv 

Mixing ratio 

141.0 

177.6 

104.3 

HCHO 

10*‘ Kg/s 

Moist convect. 

41,7 

70.1 

•28.3 

HCHO 

10*' Kg/s 

Wet dep.. 

-446.8 

-316.4 

-130.4 

HCHO 

I0*‘ Kg/s 

Dry dep. 

-4.106.2 

-2,792.2 

-1,314.0 

HCHO 

10 ' Kg/s 

Chemistry 

5,434.0 

3,555.0 

1,880.0 

HO 2 NO 2 

10 *'^ ppv 

Mixing ratio 

28.4 

35.5 

21.3 

HO 2 NO 2 

10" Kg/s 

Chemistry 

- 2,2 

3.9 

- 6.2 

CO 

10 *" ppv 

Mixing ratio 

8.3 

10.2 

6.4 

CO 

10^ Kg/s 

Chemistry 

-355.2 

-249.4 

-105.8 

H 20 

10 *" ppv 

Mixing ratio 

7,839.0 

4,272.8 

3,566.2 
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Table 2 (cont.) : Overview of global, annual mean chemical tracer budget 


Run 2 (best PI) 

Units 

Process/Conc. 

Global 

NH 

SH 

0 , 

10^ Kg/s 

Dry dcp. 

- 20.0 (81%) 

- 11.2 ( 112 %) 

-8.9 (41%) 

o. 

10" Kg/s 

Chemistry 

-0.5 (-361%) 

-0.09 (1572%) 

-0.4 (-82%) 

NO, 

10 ’* ppv 

Mixing ratio 

32.1 (12%) 

33.9 (21%) 

30.4 (1%) 

NO, 

10® Kg/s 

Dry dep. 

-31.5 (433%) 

-21.2 (589%) 

-10.2 (115%) 

NO, 

10'* Kg/s 

Chemistry 

-2,892.6 

-1,804.4 

-1,088.2 




(262%) 

(367%) 

( 88 %) 

NO, 

10*' Kg/s 

Lightning 

1,130.2 ( 8 %) 

657.1 (12%) 

473.1 (3%) 

NO, 

10 * Kg/s 

Surftce 

2,115,8 (407%) 

1,384.5 (549%) 

731.2 (138%) 

HNO 3 

10 '” ppv 

Mixing ratio 

129.3 (20%) 

136.3 (34%) 

122.2 (5%) 

HNO 3 

10® Kg/s 

Moist convect. 

-491.3 (89%) 

-243.0 (132%) 

-248.3 (46%) 

HNO 3 

10 ® Kg/s 

Wet dep. 

-189.5 (374%) 

-123.1 (529%) 

-66.4 ( 86 %) 

HNO 3 

10® Kg/s 

Dry dep. 

•475.4 (310%) 

-330.1 (402%) 

-145.4 (102%) 

HNO 3 

10® Kg/s 

Chemistry 

1,299.7 (260%) 

810.2 (364%) 

489.5 ( 88 %) 

H 2 O 3 

10 “ ppv 

Mixing ratio 

56,7 ( 88 %) 

59.5 (101%) 

53.9 (75%) 

H 2 O 2 

10® Kg/s 

Moist convect. 

-3,556.9 (56%) 

-1,667.6 (76%) 

-1,889.3 






(39%) 

H 2 O 2 

10® Kg/s 

Wet dep. 

-563.8 (89%) 

-248.9 (133%) 

-314.9 (54%) 

H 2 O 2 

10® Kg/s 

Dry dep. 

-1,672.4 (-83%) 

-904.8 (-83%) 

-767.6 (-84%) 

H 2 O 2 

10 ® Kg/s 

Chemistry 

6971.0 (63%) 

3449.0 (75%) 

3522.0 (51%) 

N 2 O 5 

10 * Kg/s 

Chemistry 

-6.3 (-1576%) 

-2.1 (-4624%) 

-4.2 (110%) 

CH 3 OOH 

10 “ ppv 

Mixing ratio 

10.0 (149%) 

9.7 (155%) 

10.2 (146%) 

CH 3 OOH 

10® Kg/s 

Dry dep. 

-328.3 (125%) 

-163.4 (122%) 

-164.9 (127%) 

CH 3 OOH 

10® Kg/s 

Chemistry 

332.5 (128%) 

166.1 (103%) 

166.5 (152%) 

HCHO 

10 “^ ppv 

Mixing ratio 

51.6 (173%) 

58.8 (202%) 

44.5 (134%) 

HCHO 

10 * Kg/s 

Moist convect. 

-13.3 (-414%) 

11,7 (499%) 

-25.1 (13%) 

HCHO 

10 ’ Kg/s 

Wet dep. 

-117.4 (281%) 

-68.2 (364%) 

-49.2 (165%) 

HCHO 

10*’ Kg/s 

Dry dep. 

-1,293.0 

-784.5 (256%) 

-508.4 (158%) 




(218%) 



HCHO 

10*’ Kg/s 

Chemistry 

1720.0 (216%) 

1010.0 (252%) 

710.0 (165%) 

HO 2 NO 2 

10 ’^ ppv 

Mixing ratio 

16.3 (74%) 

16.9 (110%) 

15.7 (36%) 

HO 2 NO 2 

10® Kg/s 

Chemistry 

-4.6 (-52%) 

-2.4 (-263%) 

- 2.2 (182%) 

CO 

10 ‘* ppv 

Mixing ratio 

2.7 (207%) 

2.8 (264%) 

2,7 (137%) 

CO 

10^ Kg/s 

Chemistry 

-72.9 (387%) 

-37.4 (567%) 

-35.5 (198%) 

H 2 O 

10 ’^ ppv 

Mixing ratio 

6,859.1 (14%) 

3,718.3 (15%) 

3,140.8 (14%) 



Table 2 (cont.) : Overview of global, annual mean chemical tracer budget 


Run 3 (PI 

Units 

Process/Conc. 

emissions) 



o. 

10^ Kg/s 

Dry dep. 

o^ 

10^ Kg/s 

Chemistry 

NO, 

lO’^' ppv 

Mixing ratio 

NO, 

10** Kg/s 

Dry dep. 

NO, 

10® Kg/s 

Chemistry 

NO, 

10® Kg/s 

Lightning 

HNO 3 

10 ” ppv 

Mixing ratio 

HNO 3 

10® Kg/s 

Moist convect. 

HNO 3 

10® Kg/s 

Wet dep. 

HNO 3 

10® Kg/s 

Dry dep. 

HNO 3 

10® Kg/s 

Chemistry 

H 2 O 2 

10 ” ppv 

Mixing ratio 

H 2 O 2 

10® Kg/s 

Moist convect. 

H 2 O 2 

10® Kg/s 

Wet dep. 

H 2 O 2 

10® Kg/s 

Dry dep. 

H 2 O 2 

10® Kg/s 

Chemistry 

N 2 O 5 

10 ' Kg/s 

Chemistry 

CH 3 OOH 

10 ” ppv 

Mixing ratio 

CH 3 OOH 

10“ Kg/s 

Dry dep. 

CH 3 OOH 

10® Kg/s 

Chemistry 

HCHO 

10 ”^ ppv 

Mixing ratio 

HCHO 

10 * Kg/s 

Moist convect. 

HCHO 

10 * Kg/s 

Wet dep. 

HCHO 

!0‘‘ Kg/s 

Dry dep. 

HCHO 

10 * Kg/s 

Chemistry 

HO 2 NO 2 

10 -*^ ppv 

Mixing ratio 

HO 2 NO 2 

10® Kg/s 

Chemistry 

CO 

10 '** ppv 

Mixing ratio 

CO 

10^ Kg/s 

Chemistry 


Global 

NH 

SH 

-19.7 (83%) 

-1 1.2 (116%) 

-8.7 (43%) 

-6.4 (-291%) 

-1.6 (-925%) 

-4.9 (-86%) 

32.0 (12%) 

33.8 (21%) 

30.3 (2%) 

-31.8 (428%) 

-21.4 (583%) 

-'10.4 (111%) 

-299.4 (250%) 

-187.8 (347%) 

-111.6 (83%) 

121.9 (0.5%) 

73.3 (0.5%) 

48.7 (0.5%) 

130,1 (20%) 

137.0 (33%) 

123.3 (5%) 

-498.1 (86%) 

-256.6 (120%) 

-241.5 (50%) 

-212.7 (322%) 

-140.8 (450%) 

-71.9 (72%) 

-496.4 (293%) 

-339.6 (387%) 

-156.7 (88%) 

1,345.6 (248%) 

842.3 (346%) 

503.4 (83%) 

65.8 (62%) 

68.7 (74%) 

63 .0 (50%) 

-3,876.5 (43%) 

-1,931.0 

-1.945.6 


(52%) 

(35%) 

-695.4 (53%) 

-348.9 (66%) 

-346.5 (40%) 

-1,855.2 (-85%) 

-967.4 (-84%) 

-887.7 (-86%) 

7748.0 (47%) 

3831.0 (58%) 

3917.0 (36%) 

-8.0 (-1263%) 

. 3.5 (-2814) 

-4.6 (-57%) 

10.9 (128%) 

10,7 (131%) 

11.2 (124%) 

-348.6 (111%) 

-170.1 (113%) 

-178.5 (110%) 

353.9 (114%) 

154,9 (118%) 

199.0 (111%) 

55,0 (156%) 

62.3 (185%) 

47.7 (1 19%) 

-2.6 (-1704%) 

11.3 (520%) 

-13.9 (104%) 

-134.6 (232%) 

-80.8 (292%) 

-53.8 (142%) 

-1,342.6 

-812.7 (244%) 

-529.8 (148%) 

(206%) 



1805.0 (201%) 

1033.0 

772.0 (144%) 


(244%) 


16.2 (75%) 

16.8 (111%) 

15.6 (37%) 

-5.2 (-58%) 

-1.2 (-425%) 

-4.1 (51%) 

2.8 (196%) 

2.8 (264%) 

2.7 (137%) 

-75.2 (372%) 

-38.7 (544%) 

-36.6 (189%) 


H20 


10*^ ppv Mixing ratio 7,804.2 (0.4%) 


4,245.7 

( 0 . 6 %) 


3,558.5 

( 0 . 2 %) 


Table 2 (cont.) : Overview of global, annual mean chemical tracer budget 


Run 4 (test sulfate) 

Units 

Process/Conc. 

Global 

NH 

SH 

0 , 

10 ’ Kg/i 

Dry dep. 

-19.7 (83%) 

-11.0 (116%) 

-8.8 (42%) 

0, 

10^ Kg/s 

Chemistry 

-6.4 (-292%) 

-1.5 (-938%) 

-4.8 (-86%) 

NO, 

10 " ppv 

Mixing ratio 

32.0 (12%) 

33.8 (21%) 

30.3 (2%) 

NO, 

10 “ Kg/s 

Dry dep. 

-31.9 (427%) 

-21.4 (583%) 

-10.5 (109%) 

NO, 

10® Kg/s 

Chemistry 

-299.4 (250%) 

-187.6 (349%) 

-111.8 (83%) 

NO, 

10® Kg/s 

Lightning 

122.1 (0.4%) 

73.3 (0.4%) 

48.7 (0.3%) 

HNO3 

10'" ppv 

Mixing ratio 

130.1 (20%) 

136.9 (33%) 

123.2 (5%) 

HNO3 

10 " Kg/s 

Moist eonvccl. 

-496.4 (87%) 

-254.9 (121%) 

-241.5 (51%) 

HNO3 

10“ Kg/s 

Wet dep. 

-213.1 (322%) 

-141.5 (447%) 

-71.6 (73%) 

HNO3 

lO® Kg/s 

Dry dep. 

-496.9 (292%) 

-340.0 (387%) 

-156.9 (88%) 

HNO3 

10® Kg/s 

Chemistry 

1,346.0 

841.6 (346%) 

504.4 (83%) 




(248%) 



H2O2 

10 " ppv 

Mixing ratio 

65.7 (63%) 

68.6 (74%) 

62.8 (50%) 

H2O2 

10® Kg/s 

Moist convect. 

-3,856.8 (44%) 

-1,924.5 

-1,932.3 





(53%) 

(35%) 

H2O3 

10® Kg/s 

Wet dep. 

-689.8 (54%) 

-343.2 (69%) 

-346.1 (40%) 

H2O2 

10® Kg/s 

Dry dep. 

-1853.2 (-85%) 

-967,1 (-84%) 

-886.1 (-86%) 

H2O2 

10® Kg/s 

Chemistry 

7712.0 (48%) 

3814.0 (58%) 

3898.0 (37%) 

N2O5 

0 

V) 

Chemistry 

-6,6 (-1509%) 

-2.3 (-4230%) 

-4.3 (-53%) 

CH3OOH 

10 " ppv 

Mixing ratio 

10.9 (128%) 

10.7 (131%) 

11.2 (124%) 

CH3OOH 

10 “ Kg/s 

Dry dep. 

-348.5 (112%) 

-170.0 (114%) 

-178.5 (110%) 

CH3OOH 

10® Kg/s 

Chemistry 

353.9 (114%) 

154.9 (118%) 

199.1 (111%) 

HCHO 

10’’^ ppv 

Mixing ratio 

55.1 (156%) 

62.4 (185%) 

47.7 (119%) 

HCHO 

10-* Kg/s 

Moist convect 

-3.4 (-1326%) 

11.4 (515%) 

-14.9 (90%) 

HCHO 

10'^ Kg/s 

Wet dep. 

-133.3 (235%) 

-80.0 (296%) 

-53.2 (145%) 

HCHO 

10-' Kg/s 

Dry dep. 

-1344,8 

-814.6 (243%) 

-530.2 (148%) 




(205%) 



HCHO 

10 * Kg/s 

Chemistry 

1806.0 (201%) 

1037.0 

•769.0 (144%) 


( 243 %) 



H 02 N 02 

10'’^ ppv 

Mixing ratio 

16.2 (75%) 

16.7 (113%) 

15.7 (36%) 

H 02 N 02 

10'’ Kg/s 

Chemistry 

-5.7 (-61%) 

-1.6 (-344%) 

-4.1 (51%) 

CO 

10 * ppv 

Mixing ratio 

2.7 (207%) 

2.8 (264%) 

2.7 (137%) 

CO 

10^ Kg/s 

Chemistry 

-74,1 (379%) 

-37.5 (565%) 

-36.6 (189%) 

H 20 

10** ppv 

Mixing ratio 

7,790.4 (0,6%) 

4,235.3 

3,555.1 


(0.9%) (0.3%) 


Table 2 (cont.) : Overview of global, annual mean chemical tracer budget 


Run 5 (test ocean) 
Ox 
Ox 
NO* 

NO* 

NO* 

NOx 

HNO3 

HNO3 

HNO 3 

HNO 3 

HNO 3 

H 2 O 2 

H 2 O 2 


H2O2 

H2O2 

H2O2 

N 2 O 5 

CH 3 OOH 

CH3OOH 

CH3OOH 

HCHO 

HCHO 

HCHO 


Units 

Process/Conc. 

Global 

NH 

SH 

10^ Kg/s 

Dry dcp. 

-20.0 (80%) 

-11.1 (114%) 

-8.9 (40%) 

10^ Kg/s 

Chemistry 

-5.2 (-338%) 

-1.0 (-1379%) 

-4.1 (-84%) 

10*^’ ppv 

Mixing ratio 

32.1 (12%) 

33.8 (21%) 

30.4 (1%) 

10® Kg/s 

Dry dcp. 

-31.6 (432%) 

-21.2 (589%) 

-10.3 (113%) 

10® Kg/s 

Chemistry 

-290.9 (260%) 

-181.4 (364%) 

-109.6 (87%) 

10® Kg/s 

Lightning 

114.7 (7%) 

67.2 (10%) 

47.6 (3%) 

10**’ ppv 

Mixing ratio 

129.6 (20%) 

136.5 (33%) 

122.6 (5%) 

10® Kg/s 

Moist convect. 

-492.2 (89%) 

-242.2 (133%) 

-250.0 (45%) 

10® Kg/s 

Wet dcp. 

-192.2 (367%) 

-124.1 (524%) 

-68.2 (81%) 

10® Kg/s 

Dry dcp. 

-479.2 (307%) 

-331.4 (400%) 

-147.8 (99%) 

10® Kg/s 

Chemistry 

1,307.4 (258%) 

814.4 (361%) 

493.0 (87%) 

I0*“ ppv 

Mixing ratio 

57.7 (85%) 

60.6 (97%) 

54.9 (72%) 

10® Kg/s 

Moist convect. 

-3,526.8 (58%) 

-1,653.5 

-1,872.3 




(78%) 

(40%) 

10® Kg/s 

Wet dcp. 

-569.4 (87%) 

-250.3 (132%) 

-319.2 (52%) 

10® Kg/s 

Dry dcp. 

-1,695.9 (-84%) 

-914.3 (-83%) 

-781.7 (-84%) 

10“ Kg/s 

Chemistry 

6985.0 (63%) 

3451.0 (75%) 

3534.0 (51%) 

10*’ Kg/s 

Chemistry 

-6.4 (-1553%) 

-2.2 (-4418%) 

-4.2 (-52%) 

l0-"ppv 

Mixing ratio 

lO.I (147%) 

9.9 (150%) 

10.3 (144%) 

10® Kg/s 

Dry dep. 

-330.7 (123%) 

-164.5 (121%) 

-166.2 (125%) 

10® Kg/s 

Chemistry 

335.5 (126%) 

164.7 (105%) 

170.7 (146%) 

10*’^ ppv 

Mixing ratio 

52.9 (167%) 

59.4 (199%) 

45.0 (132%) 

10*‘ Kg/s 

Moist convect. 

-7.9 (-628%) 

14.7 (377%) 

-22.9 (24%) 

10‘* Kg/s 

Wet dq). 

-118.0 (279%) 

-68.4 (363%) 

-49.7 (162%) 

10 ’ Kg/s 

Dry dep. 

-1,301.4 

-789.7 (254%) 

-511.7 (157%) 


HCHO 



35 


( 21 6 K) 


HCHO 

10 * Kg/s 

Chemistry 

1728.0 (2H%) 

1009.0 

719.0 (161%) 





(252%) 


HO 2 NO 2 

10**^ ppv 

Mixing ratio 

16.4 (73%) 

16.8 (111%) 

15.9 (34%) 

HO 2 NO 2 

10 " Kg/s 

Chemistry 

-4.8 (-54%) 

-2.3 (-270%) 

-2.5 (148%) 

CO 

10’* ppv 

Mixing ratio 

2.7 (207%) 

2.8 (264%) 

2.7 (137%) 

CO 

10^ Kg/f 

Chemistry 

-73.2 (385%) 

-37.3 (569%) 

. 35.9 (195%) 

H 2 O 

I0-* ppv 

Mixing ratio 

6,914.6 (13%) 

3,758.0 (14%) 

3,156.6 (13%) 
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Table 4: Downward mass transport (Kg/s) for and NO^ between lay ers 7 and 8 (approximately 1 lOmb) 

Run O, Global O, NH O, SH NO. Glnhal no mh mo c 






Table 5: Modeled OH (1x10^ molecules/cm^) and 0,(=0¥+0'D+0 j) (Tg). 
% indicates change relative to control 


00 













One standard deviation (ppb) ozone July 
















































OH (1.E4 molec/cc) run 5 (test ocean) - run 4 (test sulfate) 






















